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We adopt the general formalism, which was developed in Paper I (arXiv:0708.1233) to analyze 
the evolution of a quantized time-dependent oscillator, to address several questions in the context 
of quantum field theory in time dependent external backgrounds. In particular, we study the 
question of emergence of classicality in terms of the phase space evolution and its relation to particle 

production, and clarify some conceptual issues. Wc consider a quantized scalar field evolving in a 
constant electric field and in FRW spacetimes which illustrate the two extreme cases of late time 
adiabatic and highly uon-adiabatic evolution. Using the time-dependent generalizations of various 
quantities like particle number density, effective Lagrangian etc. introduced in Paper I, we contrast 
the evolution in these two limits bringing out key difi^erences between the Schwinger effect and 
evolution in the dc Sitter background. Further, our examples suggest that the notion of classicality 
is multifacotod and any one single criterion may not have universal applicability. For example, the 
peaking of the phase space Wigncr distribution on the classical trajectory alone does not imply 
transition to classical behavior. An analysis of the behavior of the classicality parameter, which was 
introduced in Paper I, leads to the conclusion that strong particle production is necessary for the 
quantum state to become highly correlated in phase space at late times. 

PACS numbers: 



I. INTRODUCTION AND MOTIVATION 

The subject of quantum field theory in external classical backgrounds has received significant attention in the 
literature [1-6]. Though the mathematical treatment is relatively straightforward, several conceptual issues remain 
to be settled. This paper is a sequel to a previous one in which we attempted to clarify some of these open issues in 
the context of a single oscillator with general time dependent parameters. Here, we shall shift our attention to the 
broader context of quantum fields evolving in certain time dependent external backgrounds, and in particular try to 
understand the relation of particle production with the approach to classicality in terms of phase space correlations. 

Wc begin by recalling the basic issues that were raised (and elaborated on) in Paper I [7], that one encounters in 
any semiclassical analysis: (1) In a general time dependent background, there is no well-defined procedure for defining 
the notion of particles. If adiabatic regions cannot be defined asymptotically, the usual definition of particles in terms 
of in and out states is rendered invalid. (2) A time-dependent particle concept is also required for a proper treatment 
of the question of back reaction of the quantum system on the classical background, which is not an asymptotic 
notion, many such definitions for 'particles' can be given; our aim is not to come up with a unique definition but 
a useful definition and explore its consequences. (3) The standard treatment for computing backreaction in the 
semiclassical approximation considers only the real part of the effective Lagrangian in determining the modified 
background equations of motion. The non-zero imaginary part (which directly quantifies the asymptotic particle 
content) is usually ignored in such a calculation, because the resulting equations (and solutions) will be complex (see 
for e.g. ref. [8]). This leads one to ask whether information about particle production gets encoded in some way in 
the real part. (4) The concept of particles is intimately tied with the notion of classicality of the state. This issue has 
particular relevance in studying the generation of primordial perturbations in inflationary cosmology [9-13]: in the 
standard approach, these density perturbations are assumed to be generated as vacuum fluctuations of a quantized 
(inflaton) field, but their late time evolution at super-Hubble scales is treated using classical notions, with the only 
reference to their quantum origin being in the choice of initial conditions. This naturally brings up the question of 
precisely how something that originates as a quantum mechanical entity can develop classical features over time. 

In Paper I, we confined ourselves to the case of a quantized time dependent oscillator and presented a general 
formalism to quantify the time dependent physical content of an evolving Gaussian quantum state in the Schrodinger 
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picture. In particular, wc chose to adopt an inst,ant,aneous concept of a particle, defining it in terms of instantaneous 
eigenstates specified at each moment in close analogy with a simple time-independent oscillator. (We emphasize that 
we view this as a useful definition rather than a unique definition.) As illustrated by its subsequent application to 
several toy examples, this approach provides a physically reasonable picture of the time evolution, with reference to 
particle content and phase space correlations. The main conclusions reached by our analysis of an oscillator with 
a general frequency uj{t) were: (1) Peaking of the Wigner function on the classical phase space trajectory (usually 
associated with classicality) is a rather general feature that can happen even when the oscillator evolves adiabatically, 
i.e. it is quite independent of the particle production. (2) To make the idea of interpreting classicality less ambiguous, 
we suggested considering an additional correlation measure (calling it the classicality parameter S); this quantity 
was found to track the particle creation more closely, and it does not always follow the behavior of the Wigner 
distribution. Based on an approximate asymptotic analysis of a general ui{t), we identified two distinct late time 
limits: whenever the oscillator evolves adiabatically (|e(i)| <C 1) at late times, particle creation is suppressed and the 
classicality parameter maintains an oscillatory profile; in contrast, when adiabaticity is violated (with \e{t)\ oo), 
strong particle production occurs and this is accompanied by a sharp growth in <S, whose magnitude saturates at the 
maximum value of unity. 

Here, we shall exploit the above ideas in the semiclassical context of quantum fields in time dependent external 
classical backgrounds, laying particular emphasis on making precise the connection of particle production with the 
emergence of classicality in the phase space picture. We will take up two very conventional examples, that of a 
massless minimally coupled scalar field in an FRW spacetime and a complex scalar field in a constant electric field 
background, which will prove to be quite illuminating for our purposes. (These are, of course, very well-studied 
examples which allows us one to make comparison with previous results.) The form of the action functional in these 
contexts allows, through a Fourier decomposition, to decompose the field as a set of uncoupled harmonic oscillators, 
with the time dependence of the background variables getting incorporated in the oscillator parameters. This reduces 
the field theory problem to one studied in Paper I and our formalism becomes directly applicable. (There is, of course, 
the standard issue of divergences which arise when one adds up infinite number of oscillators. While this is important 
in a different conceptual context, it is irrelevant for our purposes when we will be dealing with a quadratic action in 
a given background. Alternatively, one can imagine that a cut-off is imposed in the Fourier space at high momenta 
and attention is confined to modes well below the cut-off.) 

This paper is structured as follows. Section II provides a succinct outline of the key ideas of the formalism that 
was developed in detail in Paper I [7] to analyze a general time-dependent oscillator, placing it now in the broader 
context of the field picture. Section III is devoted to a detailed analysis, in Fourier space, of the evolution of a 
minimally coupled real scalar field in expanding FRW spacetimes and of a complex scalar field in a constant electric 
field background, with particular attention given to particle production and the phase space evolution of the Wigner 
function of the quantized field modes. Our results essentially corroborate and strengthen the conclusions that were 
drawn on the basis of the toy models we analyzed in Paper I, in relation to the dependence of phase space correlations 
on the particle production (i.e. adiabatic vs. non-adiabatic evolution). This is followed by a discussion and conclusion 
in section IV. In what follows, we shall set ?i = c = 1. 



As mentioned in the introduction, we consider here scenarios involving a quantized field, say <I>(x,t), coupled to an 
external classical background where the action for the field is expressible, by means of a Fourier transform, in the 
following form: 



Here, 5k is a real variable associated with the Fourier mode k and stands for the coordinate of the kth oscillator. 

t is the time coordinate and the dot denotes differentiation w.r.t. t. The oscillator variables mk and Wk of each 
Fourier mode incorporate the time-dependence of the external background. (To be more precise, if the background 
is characterized by time dependent parameters collectively denoted by C{t), the mass and frequency, being assumed 
to be functions of C, are implicitly time-dependent.) The form of Eq. (1) allows for a mode-by-mode analysis of the 
quantized field to be carried out, once the functions C{t) are specified. 

Although we will not go into the details here (for which the reader is referred to Paper I [7]), we provide a quick 
summary of the Schrodinger picture formalism to study the time evolution of the time dependent oscillator, to set up 
the notation. We will consider form invariant Gaussian states with vanishing mean [14-16] to describe each oscillator. 



II. FORMALISM 
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which have the general form 



^{qk,t) = A^exp[-i?k(t)«k] = ^exp 



l + 2k 



(2) 



where the wave function has been written equivalently in terms of the either the function i?k(^) or zi^(t). Using the 
time dependent Schrodinger equation, these functions can be shown to satisfy first order (but nonUnear) differential 
equations: 

2i?^ 1 

-Rk = — - - ;7mkWk (3) 
mk 2 

and 

z^ + 2iiOi,Zi,+ l(— + —){zl-l) = 0, (4) 
2 V^k mk/ 

and we shall denote the quantity ((rrik/mk) + ('^'k/'^k))/'^k- called the adiabaticity parameter, as ek(t). The problem 
of determining the quantum evolution thus boils down to directly solving for i?k or 2k (with an appropriate choice of 
initial conditions). Alternatively, since Eq. (3) and Eq. (4) are of the generalized Riccati type, they can be transformed 
into second order linear equations; for example, setting i?k = — (imk/S) {fik/ jJ-is) where Hkit) is a new function in 
Eq. (3) implies that /Xk satisfies the following differential equation: 

Mk H Mk + WkMk = (5) 

TOk 

which is same as the classical equation of motion satisfied by the oscillator variable gk, and the Wronskian condition 
AkMk — AkMk = —iW{k)/m\^, where W{\<i) is independent of time, also holds. In some cases, it turns out to be more 
convenient to solve Eq. (5) rather than start from Eq. (4). But this obviously makes no real difference as is trivially 
related to ^,k^. 

^ /oM^k + zAk V 

V^k/Uk - «Mk/ 

Although Mk(*) (or ^k) completely determines the state of the system, it is useful, particularly from the point of 
view of understanding the physical content of the evolving quantum state, to define a set of variables with physically 
reasonable interpretation, that can be built out of the wave function. For this purpose, following our approach of 
Paper I [7], we consider the following quantities: 

(1) The mean particle number (rik), that specifies for us the instantaneous particle content of the quantum state 
and is defined using the concept of instantaneous eigenstates. The mean particle number as well as the mean energy 
in a given mode (defined as the expectation value of the Hamiltonian) at any moment t are specifiable in terms of z^: 



; fk=f(nk) + ^)a;k(i). (7) 



bkP ' ■ 2 



(2) The time-dependent generalization of the effective Lagrangian Le//(k, [3, 4, 10], that specifies the transition 
amplitude between instantaneous vacuum states defined at two different moments, which again is related in a simple 
manner to Zk_: 

4 V^k mk/ 

This is to be understood as the contribution of the kth mode to the effective Lagrangian for the field. (This expression 
is somewhat formal since the functional dependence of Lgff{\s.,t) on the brackground variables is not explicit in this; 
in spite of this, we will see that it is useful.) 

(3) The Wigner distribution function defined in the q^-pk phase space of the oscillator, that provides a reasonable 
means of quantifying the approach to classicality of the quantum state in terms of its phase space evolution [10, 14, 
15, 17]. The Wigner function for the Gaussian state in Eq. (2) can be written as 



y^{qk,Pk,t) = -exp 

TT 



jL-^kW(Pk-Jk(i)gk)' 



(9) 
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where 



2^ |l + ^kP ^2|m^P 

and 

_ 2mkC^k/w(zk) _ /Wk\ 1 c^l/Xkl^ , , 

|l + 2kP " V 2 J iMkP rft ■ ^ ' 

In addition, we define an additional measure of phase space correlations, calling it the classicality parameter S, which 
is given by 



5k = ^^d=. (12) 



The classicality parameter vanishes for a vacuum state, and its value is confined to the interval [—1,1]. (The quantity 
j7k<J^ has the simple interpretation of being just twice the phase space average of (j'kPk-) O^ir analysis in Paper I 
suggested contrasting behavior of the Wigner function and the classicality parameter in different circumstances; the 
classicality parameter (but not the Wigner function) was shown to closely mirror the variation of (n). 

Since one is dealing with a field here, one can go a step further here and choose to look at the evolution in coordinate 
space which we could not do in the case of single oscillator studied in Paper I. Consider, in particular, the coordinate 
space representation of the function Zk by taking a Fourier transform and defining a function Z{x^). Working in 
the Heisenberg picture, it can be shown that z\s_ is related to the transition amplitude between: (i) a vacuum state 
defined at some initial instant ti and (ii) a two particle state (with momenta k, — k) defined using the instantaneous 
the vacuum state at time t: 

= (0,t|i-k(t)ik(^)|0,t.) 
^' (0,t|0,ii) ^ ' 

where p'k(t) = i^k(i) and A\^[t) is the annihilation operator defined at the instant t (see Paper I [7] for a derivation of 
this relation). Multiplying both sides by exp(ik • r) and integrating over k, we get the real space function: 

-^ke''-'-— 3. (14) 

This quantity has an interesting physical interpretation. If we define the positive frequency part of the field by 

A^{t)e^^---'^-—-, (15) 

(27r) 

then the state |x, y; t) = F\tl, t)F\y, t) |0, t) can be thought of as a state with two particles at x and y. The quantity 
Z{y, t) can be expressed as 

„ ,^ _ (0,^|F(x,t)f (y,t)|0,t,) _ (x,y;f|0,t,) 

This quantity can, therefore, be interpreted as the amplitude for the initial vacuum state to appear as a state with two 
particles at points x and y at the moment t. Obviously it contains some information relating to the quantum versus 
classical nature of the field. One would have expected the particles to be uncorrelated — so that Z(x — y, t) 5(x— y) 
— in the classical context, while the two particles created at two different points to retain a long range correlation in 
the quantum regime. It is therefore worthwhile to analyze the evolution of this quantity in different contexts. 

Since all the quantities defined above are related to the real and imaginary parts of Zk in a simple manner, it is 
operationally convenient to determine z\s_ and then determine the various physical quantities. Our analysis in Paper 
I thus focused on solving for the time-evolution of z in different scenarios involving various functional choices for 
ijj{t). In particular, we considered in detail toy models to illustrate the two kinds of late time evolution possible: for 
example, a frequency function uji{t) = A\/l + X^t'^, goes over into an adiabatic phase at late times {Xt » 1) with 
the adiabaticity parameter e{t) going to zero, while another function ujs{t) = A/(l + A^t^), is characterized by a 
strongly non-adiabatic late time regime (with the magnitude of e{t) diverging as t ^ co). We found that the late time 
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behavior (starting from an instantaneous vacuum state at a sufficiently early moment), while differing appreciably in 
general, shows similarities in certain respects. Although the instantaneous mean particle number follows the expected 
variation, saturating in the adiabatic example but diverging in the non-adiabatic case, as one would have expected, 
the Wigner distribution ends up peaking on the corresponding classical phase space trajectory in both situations, 
being clearly independent of the nature of the evolution, and is therefore of little use in making a statement about 
the approach to classicality. This prompted us to also look at the classicality parameter as an additional measure to 
characterize the phase space correlations; this quantity was found to closely track the particle production, and like the 
paxticle number, showed a clear difference in the two limits, remaining bounded (with oscillations) in the adiabatic 
case but saturating at unity (corresponding to a diverging q-p phase space average) in the other extreme limit. Taken 
together, these features point to potential subtleties involved in interpreting classicality in terms of the phase space 
picture. 

With a view to exploring these ideas further, we shall now proceed to take up some standard examples appearing 
frequently in field theory in the semiclassical context, which arc expected to further illustrate the two kinds of late 
time evolution mentioned in the previous paragraph. First, we will consider the often studied Schwingcr effect in 
a constant electric background, using a complex scalar field. This scenario can be brought into the form of time- 
dependent oscillator by a particular choice of gauge. When expressed in the time-dependent gauge, t his scena rio 
becomes mathematically identical to the toy model corresponding to the frequency function uji{t) = A\/l + A^t^ we 
analyzed in Paper I, which was characterized by adiabatic evolution at very early as well as very late times. On 
the other hand, the second set of examples we consider — evolution of a real scalar field in inflationary as well as 
decelerating FRW backgrounds will differ in certain respects from the electric field case, and these will be pointed 
out in the course of our analysis. It may be mentioned that the examples referred to above are very standard, and 
in which particle creation has been extensively analyzed in the literature [1, 6, 10, 14, 18]. But here we will be 
re-examining them in the broader context using the definitions introduced in Paper I which seems to lead to a more 
unified and clearer picture. 



III. ANALYSIS OF EXAMPLES 



A. Quantum evolution in a classical electric field background 



It has been known for quite a while that a classical electromagnetic field can have non-trivial quantum effects, like 
pair creation out of the vacuum (the so-called Schwinger effect) [3-5, 10]. A simple model that suffices to illustrate 
this situation involves a complex scalar field <p evolving in a spatially uniform electric background. The electric field 
E(t) can be described in a time-dependent gauge corresponding to = Az{t), with the other components being zero, 
and the action for a scalar field can then be written as [3-5]: 

A<l>] = \jd^^ [9^ + iqA^,)ct> id>' - iqAnr (17) 



= \jd^\ij dt (|(Zk|' - {kl + + {k. - qA{t)f)\q^\') (18) 

with fc^ = fc^ + fc^. Since the scalar field is complex, each mode k is associated with two degrees of freedom 
(corresponding to the real and imaginary parts of q\^); hence the number of oscillators is twice that in the case of a 
real field. 

Comparing Eq. (18) with the general action in Eq. (1), it follows that the field can be regarded as a collection of 
uncoupled oscillators, each with unit mass and a time-dependent frequency Wk(t) — -|- m? + {k^ — qA{t))^. The 
classical equation for the variable /Xk, Eq. (5), takes the form 

/tk + WkWMk = 0. (19) 

We shall choose the case of a constant electric field E > which can be described by the vector potential A{t) = 

—Et (— oc' < t < oc). This example has been widely studied as a particularly simple demonstration of particle creation 
effects in semiclassical theory, but here we are also interested in considering the phase space evolution of the quantized 
field. 

By the transformation r = kzl^~qE -\- ^JqEt, Eq. (19) can be converted to a simpler form: 

^ + (r^ + Ak)/Xk = 0, (20) 
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where Ak = (fcj^ + m^)/qE. This corresponds to an oscillator with the time-dependent frequency ciJk(T) = + Ak- 
For this frequency function, the adiabaticity parameter given by ek(''') = t /[\\^ + r^)^/^ vanishes in the asymptotic 
limits (r — > ±00) , so in and out vacuum states can be defined in the adiabatic regions. If one starts out in the in 
vacuum state (defined by (rik) ^ as r ^ —00), the non-equivalence of the two asymptotically defined vacua implies 
that the state will appear populated with quanta measured with respect to the out vacuum at late times. 

As a first step, we would like to understand how the instantaneous particle number, introduced in Paper I, evolves 
with time. This quantity is, of course, expected to coincide with the definition based on in-out states in the adiabatic 
regions. A computation of the particle content requires one to solve for 2;k or /Zk and, for the constant electric field 
case, Eq. (20) is identical in form to the Schrodinger equation for an inverted harmonic oscillator, which is exactly 
solvable. We can thus obtain the analytic solution, and this will allow us to determine all the physical properties, 
including the exact asymptotic particle content at late times. 

The linearly independent solutions of Eq. (20) are the parabolic cylinder functions iD,y»(±(l + i)r) and their complex 
conjugates [19], where fk = —1/2 + iAk/2. The particular solution depends on the initial condition chosen. If it is 
assumed that the field modes start out in the instantaneous vacuum state (vanishing particle number) in the asymptotic 
past (r — > —00), then the wave function of each oscillator must satisfy 



.lim. Vk(r) ^ (''^)'''exp (-"^g^ - ^ / ..dr ) (21) 



which translates into the following condition on ^,\^■. 



lim ^ « iu;k(T) (22) 

where the prime denotes differentiation w.r.t. r. {qE > has been assumed, so t and r have the same sign in the 
asymptotic limits.) The particular solution of Eq. (20) which has the required behavior is [— (1 + i)T]. This may 
be seen by using the asymptotic form of this function in the r — > —00 limit: 



lim D,.[-(l + i)r] « {V2\T\)<e^-'^-"^ (23) 

r— ^ — 00 ^ 

which gives 



iMr)-7;^ (24) 



^k |t| 2a;k(T) 

and the Wk/wk term vanishes in the r —00 limit. We next need to obtain the behavior at late times, for which the 
r — > 00 limit of /ik is required. This is given by 

£',.[-(l + i)r] ™ e-''^'^i(y2T)^^eT-i:-*4 ("l + i^iM^ 



e 2 1. 



''^+i)(V2r)-e-^l-+^^ {l-i^"*^^^lX^^^^) 



where 

A^ = {^f2)<e-'"-i< B^ = -i-^l^-{y/2Y-e-''i(^<+--^ e^(^) = ^ + ^ ln(r) (26) 

i fkJ ^ ^ 

4 4 
Using this approximation, it can be shown that in the late time limit, 

^k « -^j^^ " 4^ + ^^k) + ^^e^'^--^^)) - 27^ke^^^^(^) (c, - (2 - n^e^^<>^^^^))) + O {r'^) 

(28) 



and 
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where TZu = Bj^/A^. Using the above asymptotic form for z^, the late time limit of (rik) can be determined: 

lim (nk) = lim ^'f = (29) 



Substituting the exact expressions for ^k and Bk gives 
and hence 



= S = = 2e-^^/^cosh ( ^) (30) 



(^2 _j_ ^2 \ 
-TT " ^ . (31) 

qrj ) 

This expression is, of course, identical to the standard result which is normally computed using in-out states [3, 5]. 
This agreement in the asymptotic limit shows that our choice of the time-dependent particle number to be a reasonable 
one. 

Given the exact solution, one can obtain the mean particle number as a function of time in terms of parabolic 
cylinder functions. The analytic forms are not illuminating but the plot giving the variation of the mean particle 
number as a function of the dimensionless variable r — given in figure 1(a) for the choice of Ak = 1 — has several 
interesting features. The particle number (nk) is nearly constant at very early and very late times as expected. It 
starts from zero in the asymptotic past, because of our specific choice of initial state. At late times, in the r ^ oo 
limit, it saturates at the constant value given by Eq. (31). It may be noted that since the variable r actually depends 
on both the original time coordinate t and the momentum fc^, the value of r for a given t varies from mode to mode, 
implying that modes with different momenta would saturate and reach the adiabatic phase corresponding to r » 1 at 
different times. More importantly, at intermediate times (when the adiabaticity parameter e is appreciably non-zero), 
(rik) is characterized by oscillatory behavior superimposed on a steadily increasing mean, which can be interpreted 
as a continual interplay of creation and annihilation of quanta. The amplitude of these oscillations gets progressively 
diminished as the evolution proceeds into the late time adiabatic regime, and the mean settles down at a constant 
value. (Though Schwinger effect is extensively discussed in the literature, the authors are not aware of a plot like the 
one in figure 1(a) in any of the works.) 

A handle on the nature of the time- evolution of the oscillator can be gained by directly visualizing the evolution of 
Zk in the complex plane. This is shown in figure 1(b). The complex trajectory of Zk starts from the origin, and after 
meandering around for a while in the intermediate phase, ultimately ends up, for large values of r, on a |^;k| = constant 
curve, corresponding to an asymptotically constant mean particle number. 

When we study particle production in cosmological backgrounds — especially in the context of inflationary models 
— it is usual to focus attention on the power spectrum, which is related to the Fourier transform of the two-point 
correlation function for the field. To make latter comparison between Schwinger effect and cosmological particle 
production easier, we shall compute this for the Schwinger effect. The power in the kth mode is just the variance of 
^k in the given quantum state, and it can be written as 

,2. _ 1 _ (2(nk) + 1) , ((nk) + l) . . . . 

^'-^ - 2(i?k + i?£) - 2.^ + (32) 

The power, though related to the time dependent mean particle number, is not completely expressible in terms of it, 
being also dependent on the phase of Zk (the particle number however involves only the magnitude of Zk). Only under 
the condition of zu_ taking on a constant value, are the two quantities directly connected. In the present situation, z^. 
remains oscillatory in the adiabatic limit at late times while the particle number approaches a constant value, so the 
power, for large r, is approximately given by 

ill) - -^^] + Re (33) 



which can be written as 



2t J \ T J \ Tlk 



1 



(9k> « z ( e"''^" + 2 + V e-'^^" (e"'^^- + 1) cos (r^ -h Ak In r -h ru) ) (34) 



■ + rk)) 
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where 72.k = |7?.k| exp(irk). If the third term in the brackets, which is rapidly oscillating, is dropped, then the variance 
at the leading order is just 

(ql) « I ((^^k> + (35) 

and the power in the adiabatic limit can be related to the mean particle number. But this connection is not quite clear- 
cut for two reasons. First, it requires our ignoring the additional oscillatory term which allows us to ignore the phase 
information contained in z^. Second, the correlation actually dies down as r when r — > oo and its interpretation 
is unclear. This point will be taken up again when making a comparison with the field theory in the de Sitter 
background. 




FIG. 1: Plots of the particle number (rik) and the complex trajectory of with time for the constant electric field case 
(tJk(T) — V-^k + T^) for Ak = 1. (wk) saturates at a finite value at late times, Zk starts from zero, and at late times ends up 
going around in a circle centered on the origin of the complex plane. 

Let us now turn to considering the question of quantum versus classical nature in terms of the phase space evolution 
of the Wigner function given by Eq. (9), particularly in the asymptotic limits. The relevant functions which determine 
the Wigner function, and Jk are shown in figure 2. 
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FIG. 2: Variation of the fmictions \/qEa^ and J-^j \JqE for the electric field case (for Ak = 1). The Wigner function starts 
uncorrelated and peaked on the p— axis, and again at late times ends up peaking on the p— axis. 

As is clear from these plots, the Wigner function starts uncorrelated and sharply peaked on the vertical pk axis. 
At late times it again ends up being concentrated on the pk axis. This behavior may be compared with that of the 
corresponding classical phase space trajectory, which, in either asymptotic limit (|r| 3> 1) has the approximate form 



1 



(Ak + r2)iA 



rRe 



Ce 



iT y/^k+r''-'dT' 



1 



1 



4r2 



Re 



Ce 



(36) 
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where C is an arbitrary complex constant, with the momentum being 



Pc 



qc 



qc_ 
'2t 



4|r|3/2 



Re 



(37) 



It is clear that in both asymptotic limits, qc ^ and the trajectory is along the axis. In other words the Wigner 
function is concentrated on the classical trajectory both at very early times (when the field is in the in vacuum) and 
at late times (when particle number has reached an asymptotically constant limit. Hence, this aspect (i.e. peaking 
on the classical trajectory) of the evolutionary behavior of the Wigner function alone cannot be used in a direct 
manner to conclude the emergence of classicality. Obviously interpreting classicality merely in terms of peaking on 
the classical trajectory is suspect, since this can happen even when the oscillator is in a near- vacuum state (which 
happens to be the case here at early times). 

It would be of interest to compare the above behavior with that of the classicality parameter iSk defined in Eq. (12), 
for this situation; its time variation is depicted in figure 3. This correlation measure clearly displays an asymmetry in 
time and more closely tracks the particle number: 5k stays very close to zero in the course of the early time adiabatic 
evolution, but in the late time adiabatic regime, ends up oscillating about 5k = 0. The oscillatory nature of 5k 
refiects the oscillations of the imaginary part of Zk, as it circles around; this becomes obvious by recalling the relation 
connecting Zk with 5k. As is evident from the plot, the time-averaged mean of this oscillatory variation is nearly zero, 
but the variance of 5k has a non-zero finite value that stays nearly constant in the late time limit. Thus, although 
the Wigner function provides little help in addressing the question of classicality in this situation, the variation of 
5k points to the state becoming more correlated at late times, basically due to some amount of particle creation 
occurring over time, in comparison with the nearly zero value at early times. Nevertheless, none of these measures 
indicate a clear emergence of classicality in this context. We will see later that the situation is quite different in the 
case of evolution in the de Sitter background, showing a sharp difference between the Schwinger effect and evolution 
of a quantum field in a deSitter background. 




FIG. 3: Variation of the classicality parameter cSk with time for the constant electric field case (Ak = 1). 5k is nearly zero in 
the asymptotic past, but ends up oscillating about zero with a finite amplitude at late times. The variation of 5k reflects the 
complex trajectory of Z]^, which at late times is a circle centered on the origin of the complex plane. 

One can also determine the form of the time-dependent effective Lagrangian L^ff in the asymptotic limits from 
Eq. (8). Using Eq. (24), it can be shown that Zk has the following approximate form in the adiabatic early time limit: 



« --6k(r) + 0(e2) 
where ek(''") = uj')^/lu^, and the mean particle number is given by 



4(r) 
16 



The asymptotic form of Zk allows one to derive an approximation for ie// in the r 
consequently be written as 



ReL 



^k('r)e^(r) 
16 



(nk)wk(r) 



(38) 
(39) 

-oo limit. The real part can 
(40) 
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to lowest order in e. This expression has the simple interpretation of being the energy drained from the background 
electric field due to particle production, — <x'k('^k)- (Recall that the effective Lagrangian and effective Hamiltonian 
differ by a sign in this case.) Thus, to lowest non vanishing order in e, the real part of the effective Lagrangian 
incorporates information about particle creation during adiabatic evolution starting from a vacuum state, and this 
suggests that to lowest order in the adiabaticity parameter, the backreaction due to the production of particles can 
be incorporated using the first term in HcLeff. 

As for the late time adiabatic Hmit of large r, the asymptotic form of Zk in Eq. (28) suggests that at the lowest order, 
the real part of the effective action (per mode) will pick up dominant contribution over the late time adiabatic phase 
of evolution from the non-oscillatory term of magnitude in the expression for ReLg//- The imaginary part, 

on the other hand, is directly related to the asymptotic particle content with respect to the out vacuum state. One 
can further obtain an expression for the effective action for the field itself, by using the result that the contribution 
coming from one mode is just (1/4) ln(l + (rik)) (see Eq.(38) of Paper I). We get: 



= lSl„(l + e-) = |t (41) 

k •> \ / n=l 

To arrive at the last expression we have summed over every k twice, since each mode of a complex field is associated 

with two degrees of freedom, and used the Taylor series representation of ln(l + e~'^^). The summation over k has 
been transformed into an integral by the standard procedure of box normalization with periodic boundary conditions 
over a (large) volume V. The integrand has no dependence on the z-component of the momentum; so making the 
change / d^k ^ J dkz J 27rA;„dfc„ and integrating over fc„ gives 

" J qE n2 (27r)3 ^ ' 

n=\- 

Making the replacement kz/qE T (regarding T as a time coordinate) brings the above expression into the form 

W,„ = VTf;tlpi|2!e-'?f^ (43) 

which is of course the same as the standard expression obtained using in-out states [3, 5]. 

To summarize, our analysis of the time evolution in a constant electric field, which maps to a harmonic oscillator 
problem with late time adiabatic evolution, reproduces all the standard results. In addition it indicates that the 
peaking of the Wigner distribution on the classical phase trajectory to be a fairly general feature, quite independent 
of the amount particle production. It occurs even in the course of adiabatic evolution and in the initial quantum 
vacuum state. The q-p correlation, in contrast, remains bounded and oscillatory in the late time adiabatic regime 
where particle production is suppressed, and is evidently more sensitive to the nature of the evolution, viz. whether it 
is adiabatic or non-adiabatic. The relevance of these results will become more apparent when we study the evolution 
in cosmological background. 

We shall now take up another scenario, that of expanding spacetime backgrounds, and consider time-evolution 
of a quantized field in specific models which illustrate both the kinds of late time evolution: adiabatic as well as 
non-adiabatic. 



-TrnA 



B. Quantum evolution in expanding universes 



The Friedmann model provides an instructive setting for analyzing quantum effects in a time-dependent gravitational 
field. A spatially flat Friedmann universe [1, 9-12] is described by the line element 

ds^ = df - a^{t) d-x^ = a^{ri) (dr]^ - rfx^) , (44) 

where t is the cosmic time. a(t) is the scale factor and rj denotes the conformal time with dr] = dt/a{i). We will 
consider a massless and minimally coupled real scalar field $ in this spacetime (which, for example, could even describe 
linearized tensor perturbations [1, 9, 10, 12]), governed by the covariant action 

A[^] = i / d'^x^^^,^^^'^. (45) 
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The homogeneity and isotropy of the Priedmann metric (44) allows one to decompose the scalar field $ as 



(3k(??) e' 



ik-x 



(46) 



(27r)3/2 



The Qk are in general complex and can be written as Qk = cik + ibk, Ok and 6k being real. But since $ is a real 
field, Q-k = <5k t;his imposes the constraints a_k = flk and 6_k = —6k, which means there actually is only one 
independent degree of freedom for every k, and this can be chosen to be Ok for one half of fc-space and 6_k for the 
other half. Relabeling these real variables as over all of A;-space, the action can then be expressed as 



where the primes denote differentiation with respect to the conformal time 77 and k = |k|. The action (47) describes 
a collection of independent oscillators with time dependent mass m = a? and frequencies u = k. The quantized field 
can thus be analyzed mode-by-mode, and the general formalism developed to study a time dependent oscillator can 
be directly applied here. 

Let us now quickly recall some standard notions. In an expanding universe (described by a scale factor that is an 
increasing function of the cosmic time t), all physical length scales get stretched by the factor a(t). Depending on 
the equation of state of the matter source that drives the expansion, the Friedmann model can undergo acceleration 
(characterized by d{t) > 0) or deceleration {a{t) < 0). The scale that determines the dynamics of a field, say, in 
such a background is set by the Hubble radius, defined as Rnit) = H{t)~^ = {a/a)~^ = {a' /a?)~^. Any given field 
mode with a comoving wave vector k will have its physical wavelength red-shifted by the a{t) factor, and there are 
two possible regimes: one, when the proper wavelength is smaller than Rh (sub-Hubble scales) and the other when 
the mode is 'outside' the Hubble radius, its physical scale exceeding Rh (super-Hubble scales). The dynamics of the 
evolution of field modes in these two limits can be qualitatively quite different. 

The adiabaticity parameter for any given mode in this setting has the form 6(77) = 2a'/(a[k[) (or in terms of cosmic 
time, e{t) = 2a/|k| oc X{t)/ Rnit)), which is essentially the ratio of its physical wavelength to the Hubble radius. 
This immediately tells us that in an accelerating universe (characterized by a > 0), the adiabaticity parameter is an 
increasing function of time, so any given mode is expected to make a transition from initially adiabatic to highly non- 
adiabatic evolution at sufficiently large times, following Hubble exit. This trend will reverse in the case of decelerating 
models (a < 0) and any mode will eventually end up evolving adiabatically. 

These features clearly show that this scenario of a quantum field in an expanding gravitational background is quite 
different (conceptually and algebraically) from that of the constant electric field case discussed earlier. In the case 
of electric field the evolution was adiabatic both at early and late times allowing us to define asymptotic vacuum 
states. Rigorously speaking, we can define an adiabatic in-vacuum at early times in the case of accelerating models 
while one can define an adiabatic out-vacuum at late times in the case of decelerating models. In such situations, the 
time-dependent definition of the particle concept is expected to provide a physically reasonable alternative, and will 
be adopted in our analysis of the quantized field in the examples which follow. 

On general grounds, particle creation effects are likely to be significant only in the non-adiabatic, i.e. super-Hubble 
phase of evolution (where ek(??) = k/aH > 1). The expectation of strong late time particle production in inffationary 
universes, contrasted by saturation of the particle content as well as energy in the decelerating models at large times, 
naturally follows from this. (We again stress that terms like 'particle production' are used here in a specific sense, 
based on the idea of instantaneous vacuum, and should not be misinterpreted.) 

The standard approach to determine the evolution of the system starts with the equation for /Xk, Eq- (5). (Isotropy 
of the background means that all quantities will be functions of only |k| = fc, and this allows us to replace the label k 
by k; we shall stick to this notation from now on.) The solutions of Eq. (5) are usually straightforward to obtain, and 
once the particular solution for a given initial condition has been chosen, it completely fixes the quantum state. The 
table below lists the various physical quantities expressed in terms of Zk (which is directly computable from /x/j): 




(47) 
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The wave function 


V.fe,r,) = exp 


Mean particle number 




Mean energy per mode 


£fc = ((nfc) + 5) 


Wigner function W {qk jPkt't) 




Classicality parameter 


2/m(sj,) 


The power spectrum 


= ^'^.'=i+^^ + "^*2+^'««(^^) 


Effective Lagrangian 


ImLe// = i ^ Jie(;^fc) = J i| ln(l + {uk)) 



Our basic interest is in understanding how the time-dependent particle content and the Wigner function describing 
the phase space correlations of the quantum state evolve in the expanding spacetime background in various limits. 
To keep the analysis clear and simple, we shall consider cosmological models characterized by power-law expansion 
with the scale factor having the form a(t) cct^.p lying in the range [0,oo]. This set covers accelerating (p > 1) as 
well as decelerating (0 < p < 1) types. As mentioned before, our intention is to consider examples that illustrate 
the two possible extreme cases of late time evolution, and the power law models prove adequate for this purpose. 
In particular, the accelerating models have close correspondence with the toy model of the frequency function W3 
we considered earlier, which is non-adiabatic at large times, while the non-inflationary ones share the feature of an 
asymptotically adiabatic phase in common with two other toy models (those with the frequency functions ui and 
0^2). We shall begin by picking one representative example from each of these two distinct groups, the de Sitter 
model and the radiation-dominated model respectively, and try to bring out the essential features that distinguish 
(or are shared in common by) them. Following this, a model combining an initial decelerating phase with a late time 
accelerated period of expansion will be analyzed, and finally a general approximate asymptotic analysis will be carried 
out separately for inflationary and decelerating models. 



1. The de Sitter universe 



The de Sitter universe [9-12] has exponential a{t) and can describe the inflationary universe. It essentially corre- 
sponds to taking the limit of the index p ^ 00 in the power law models with suitable normalization. In terms of the 
conformal time r], the scale factor is given by 

«(^) = -^ (48) 

where rj <0. The adiabatic parameter Ck = '^/k\r]\ = 2aH/k clearly blows up as the scale factor goes to infinity. This 
naturally means that every mode, starting out at sub-Hubble scales at early times, will 'exit' the Hubble radius and 
continue to remain super-Hubble thereafter. As indicated earlier, the divergence of ekiv) indicates the absence of an 
adiabatic out- vacuum state, and hence a natural definition of particles at late times. Consequently, and definition for 
particle is intrinsically ad hoc and one need to base the discussion on something which is useful. We will see that the 
time-dependent instantaneous particle concept is a fairly reasonable choice in such a situation. 
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To determine the time evolution of the quantum state, we will begin with Eq. (5). For the expansion given by 
Eq. (48) the general solution for (j.^ can be expressed in terms of simple functions [19, 20] as follows: 



M/s(a) = 



1 



2ka 



1 



uiH 



^ + nk\i + 



iaH 



(49) 



where TZk is an arbitrary complex number. If one assumes the mode to be set in the instantaneous ground state at 
the moment a = Ui, this corresponds to imposing the condition Rk{a,i) = fcaf/2. The function TZk then turns out to 
be 



2ik 



exp 



2ik 



(50) 



Here, we shall choose to work with a state that begins in an instantaneous vacuum in the asymptotic past; this 
implies a, 0, so that TZ^ vanishes for all k. This state is generally known in the literature as the Bunch-Davies 
vacuum [1, 9, 10], and its evolution will be considered in the analysis which follows. Once the form of Hk is specified, 
it is trivial to show that 



Zkia) = 



2k 

aH 



2ik 
aH 



(51) 



which is quite simple. From this it follows that the time-dependent mean particle number is given by the following 
expression: 



{rik) 



\zk (a) 



kfe (a) 



4fc2 



Rh 



(52) 



Ap being the physical wavelength of the mode. The particle number sharply diverges in the late time non-adiabatic 
regime, i.e. at super-Hubble length scales {k/aH <^ 1) as the accelerated expansion drives rapid particle production. 
The time- variation of (n^) as well as the trajectory of zj. in the complex plane for the Bunch-Davies state are depicted 
in figures 4(a) and 4(b). Zk starts from the origin, and eventually ends up at the limiting point corresponding to 
Zk = +1- This is in sharp contrast to the evolution in the case of Schwinger effect, in which the z goes around a 
closed loop in the complex plane at late times (see figure 1(b)). 

It may be noted that since the mean particle number here grows monotonically with time, one can trade off the 
time dependence for dependence on (rife), and therefore express Zk entirely in terms of the mean number of particles: 



Zkiirik)) 



1 + 



2ik 
aH 



\/{nk) 

i + \/ {nk) 



(53) 



This in turn allows one to express the wave function, too, explicitly as a function of just the mean particle number: 



[nk) 



4)fc3 



ttH^ \1 + A{nk) 



1/4 



exp 




(rik) 



nk) 



4 



(54) 



This is possible here only because of a one-to-one relation between a and (rife). Such a transformation will not be 
unique in general (for e.g. if the mean particle number is oscillatory, the inversion would give multiple values of time 

for the same (rik)). 

In the literature, what is normally computed for infiationary models is not the particle number but the power 
spectrum of the field [9-11, 13], and we would like to connect it up with the particle concept used above. As 
mentioned earlier, the power in the kth Fourier mode is just the variance of the qk variable in the given quantum 
state, and is expressible as follows: 



ill) 



1 

2 {Rk + Rl 



(2(nfe) + l) , ((nfe) + l) 



2A;o2 



+ 



ka? 



-Re(^fe). 



(55) 



From Eq. (55), it is clear that the power spectrum is not in one-to-one correspondence with the particle content, since 
it also depends on the phase of Zk through the Re(2;fe) term. It is only under certain conditions — like when Zk goes to 
a constant value — that a direct correspondence between the two can be established. This is possible, for example, if 



14 




FIG. 4: Plots of the particle number (n^) and the complex trajectory of zt with the dimensionless variable aH/k for the 
Bunch-Davies state in the de Sitter model. Sub- and super-Hubble regimes are color-coded blue and red respectively, {nk) 
starts from zero and diverges at late times. The trajectory of z^, starting from the origin, remains confined to a quadrant in 
the complex plane, and at late times (super-Hubble scales) approaches the point zu = 1- 



the power is evaluated at super-Hubble scales {k/aH ^1), because in this limit, — ^ 1, and from this it also follows 
that the particle number is large ((n^) 3> 1). Consequently, we have 

(ql) « 2(n,)/(fca') (56) 

and thus, at large scales, the power in any given mode does provide a measure of its particle content. Comparing 
Eq. (56) with the corresponding equation Eq. (35) in the case of Schwinger effect, we find that there are significant 
differences. First, the situation is mathematically much clearer in the case de Sitter universe in which the asymptotic 
limit of z is unambiguous. This is not the situation in the case of Schwinger effect and one need to drop oscillatory 
terms in a somewhat arbitrary fashion to arrive at Eq. (35). This can be tracked down to the fact that in the case of 
de Sitter universe the late time evolution is highly non-adiabatic while in the case of Schwinger effect the late time 
evolution is adiabatic. As shown in Paper I, the adiabatic evolution generally leads to a late time oscillatory behavior 
for Zk making the relation between power spectrum and particle content ambiguous. 

For the case of de Sitter evolution starting from the Bunch-Davies state, the power is given by 

<''">-£('+^^)- 

When evaluated for every mode at Hubble exit {k/aH = 0{1)), this expression yields k^{ql) cx H^, corresponding to 
a scale-invariant spectrum. 

We now turn our attention to understanding the quantum versus classical nature of the state. We will start with 
the evolution of the Wigner function in phase space. The Wigner function is determined in terms of the two variables 
cr| and Jfc which are given by 



Figures 5(a), 5(b) show the variation of these functions with time. From the above expressions, it is clear that 
Jk — > 0, (T^ — > oo at early times {k/aH — > oo) i.e. sub-Hubble scales, which corresponds to a state that starts sharply 
concentrated on the g-axis. At late times {k/aH 0), however, J'k ^ oo with a finite a'^, which causes the Wigner 
function to end up peaking on the p-axis. 

Let us compare the above behavior with that of the corresponding classical phase space trajectory. For the case of 
exponential inflation, the general solution for qk in terms of the time coordinate rj can be written as 



qk = Rc 



C,(l + -)e-'' 



(59) 
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where C is an arbitrary complex number. Writing C = Ce*'^, this solution can be re-expressed as 

C 

qk — Crj cos ikrj + c) — sin [kij + c) . (60) 

k 

The conjugate momentum = {a^ corresponding to the above is then given by 

Pk ^ - {Ck/H^ri) sin [h] + c) . (61) 
The classical trajectory of the system in the phase space thus takes the form 

* = -^Vk ^c^[\- -^pIJ (62) 

where rj is to be understood as being an implicit function of (which can be found by inverting Eq. (60)). At late 
times (ry 0) i.e. at super-Hubble scales, {qk/Pk) 0, indicating a trajectory that ends up along the vertical pk 
axis, and this is precisely what the Wigner function gives. As for the other limit oi r] —f — oo, we get Pk with 
rjpk remaining finite, which obviously gives a trajectory along the qk axis. This again matches with the analysis 
based on the Wigner function. Thus, the Wigner function is found to peak on the classical trajectory not only in 
the non-adiabatic late time phase when strong particle creation takes place, but even at very early times when the 
quantum state of the mode is in fact assumed to be close to a vacuum. This may be taken to suggest that this 
aspect of the evolutionary behavior of the Wigner distribution is not sufficient as far as providing information about 
classicality is concerned. We reached a similar conclusion in the case of the Schwinger effect. 




FIG. 5: Variation of the dimensionless functions k^a^/H^ and H^J'k/k^ for the Bunch-Davies state. Sub- and super-Hubble 
regimes are color-coded blue and red respectively. The Wigner function starts uncorrelated and concentrated on the g-axis; at 
late times (aH/k — > oo), q/p — + representing a state concentrated on the p-axis. 

In contrast, the classicality parameter gives a more meaningful description here. It can be directly determined from 
Eq. (58), and is given by 

Sk{a) ^ _ , -2,0^. (63) 



^1 + [aH/kf Vl + 4("'=)' 



Let us look at the asymptotic limits of this expression. At sub-Hubble scales {k/aH — > oo), the mean particle number 
is small, and 



Sk « - [1 + O {{uk))] (64) 

while at late times, as (n^.) grows sharply at super-Hubble scales {k/aH 0), 

Sk « -1+^ (65) 
Sink) 



to the lowest order. (At Hubble crossing, k/aH ~ 1 and so Sk ~ l/\/2-) It is clear that this quantity Sk follows 
more closely the particle content, in that it remains close to zero in the initial adiabatic sub-Hubble phase, but grows 
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subsequently, saturating at unity in the large {uk) limit. (This, it may be remembered, was also found to happen 
in the non-adiabatic toy model corresponding to 0^3 considered in Paper I [7].) Thus, if Sk is used as a measure of 
the degree of correlation, then one can unambiguously conclude that the state is evolving towards classicality, which 
agrees with the standard lore in this subject. But note that the Wigner function peaks on the classical trajectory 
both at early and late times and hence is not such a good descriptor. 
The plot of Sk for the Bunch-Davies state is shown in figure 6. 




FIG. 6: Plot of the classicality parameter Sk for the Bunch-Davies state in the de Sitter model. Sub and super Hubble regimes 
are color-coded blue and red respectively. \Sk\ starts from zero in the asymptotic past and saturates at unity at super- Hubble 
scales. 

The contribution to the time-dependent effective Lagrangian coming from the kth mode also can be computed 
using the expression for Zk, and is given by 

* "'4 A iaH\ 

where the real part can be rewritten in terms of the mean number of particles in the mode k, as 

In the asymptotic past {k/aH 00), the imaginary part goes to zero, while the real part gives 

ReLe// « fc(nfe>+0((n,)2) (68) 

(recalling that the mean particle number in the Bunch-Davies state vanishes in this limit). To lowest order in (nfc), 
this can be interpreted as the energy drained from the gravitational field that goes into creating particles, ~{nk)k. 
This again suggests, as in the electric field case, that information about particle production during adiabatic evolution 
starting from a vacuum state (for the quantum field) is encoded in the real part of the effective Lagrangian, and that 
one can take into account backreaction due to the produced particles using ReLeff. This makes sense and lends 
further credence to our time dependent definitions. 

In the other limit of late times (corresponding to k/aH 0), the imaginary part blows as a consequence of strong 
particle creation, while the real part assumes the approximate form 

ReLeff ^k(^l- (69) 

which saturates at the limiting value ReL^ff = k asymptotically. Clearly, in this limit where (rij.) 3> 1, the corre- 
spondence of ReLeff with particle creation which holds in the adiabatic regime is not valid, and the interpretation of 
the above expression is presently unclear. 

Finally the physical content of the quantum state can also be clarified in terms of the spatial Fourier transform of 
Zk, which, as was shown earlier, can be thought of as the amplitude for a pair of particles to be produced from the 
vacuum at any given instant. This function Z(r,a) is given by 



d^k 1 r°° k sinikr) 

— / 



Z(r,a) = / e-"-z,— = / ^^^-^dk (70) 

aH 
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which can be written in terms of standard functions: 

, i (a^E'^\ r 1 \ ( ..H ( raH\ ran ^ fraH\ , 

Ei{x) being the exponential integral function [19]. While this is not very illuminating, it has clear asymptotic forms 
which confirms our interpretation. When a ^ 0, the second term in the square brackets above vanishes and 

On the other hand, when a — > oo it reduces to a Dirac delta function: 

lim Z{v,a) = 5{v). (73) 

a— ^oo 

(While one can get this result by a careful limiting procedure, it is obvious from the fact that at late times, ^ 1.) 
The first limit found above can be interpreted as suggesting that at early times, the particles created at two points are 
strongly correlated, implying an underlying quantum nature, and the correlation falls off with separation only as r~^. 
The large a limit of a delta function, on the other hand, implies particles with non-zero separation are uncorrelatcd, 
in line with the classical notion of a particle as a spatially localized entity. Thus, Z{r, a) provides a measure of the 
classicality and particle content of the state consistent with what we found earlier. 

To recapitulate, the de Sitter background provides an example in which the late time evolutions is highly non- 
adiabatic with no sensible out- vacuum in the conventional terminology. Our formalism, nevertheless produces sensible 
results consistent with what have been known in the literature and shows that the state evolves towards 
classicality. This is clear in the behavior of the classicality parameter (which increases from zero to unity in the course 
of evolution) and that of Z(r, a) which becomes a delta function at late times. The Wigner function, on the other 
hand, is not a good descriptor of classicality since it is peaked on the classical trajectory both at early and late times 
and does not track the classical trajectory in the intermediate epochs. This shows a contrasting behavior vis-a-vis 
the evolution in a constant electric field background which leads to adiabatic evolution both at early and late times 
allowing for adiabatic in-vacuum and OTit-vacuum to be defined. The evolution of the system towards classicality is 
much less clear in this case. Altogether, these results closely resemble those found for the non-adiabatic and adiabatic 
toy examples of Paper I [7]. 

We now move on and choose an example of a decelerating universe, to illustrate the other alternative of late time 
adiabatic evolution. We expect this to present a contrast with the behavior found above and resemble the Schwinger 
effect more closely. 



2. The radiation-dominated model 



We shall now work in the setting of radiation dominated expansion [11, 12], that corresponds to setting p = 1/2. 
The radiation dominated model is thus governed by the expansion law 

a{r])=aoil (0 < r? < oo). (74) 

Here, the adiabaticity parameter given by e{r]) = 2/kr] = 2ao/{ka) steadily drops to zero with expansion, and particle 

production in any given mode is naturally expected to be strongly suppressed in the late time sub-Hubble phase of 
evolution. In fact, in being characterized by a late time adiabatic regime, the radiation dominated model bears close 
resemblance to the constant electric field case that was considered earlier and contrasts with the de Sitter example. 

For the linear scale factor in Eq. (74), a" = and the solutions of the classical mode equation Eq. (5) are particularly 
simple, being just plane waves; the general solution for /Xfc can be written as 

Mfe(a) = ^ (e-'^o + TZke-^'-^o) (75) 



2fco 

7?.fe being an arbitrary constant coefficient. If it is assumed that the mode starts off in the instantaneous vacuum state 
at the moment a = a,, then 7?./; = — (1 + 2ikai/ao)~'^ ^xkaijaa ^ o^^^j jg straightforward to determine the form of Z}^: 

2ifc(^)7efce-^''=4 + (l + 7e,e-^^'=4) 
^^(«) = ( \ ( ■ (76) 
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If the state is extended into the past to a < a^, ^ — 1 and the particle number is found to diverge as the initial 
singularity is approached, for arbitrary TZk, and it is not possible to define an instantaneous vacuum state at = 
when the spacetime is singular. Therefore, we have to confine ourselves to states which start off in the ground state at 
some Qi > 0. The behavior of the field mode prior to will not be of relevance to us, and we will limit ourselves only 
to the subsequent evolution of the wave function (for a > ap). This situation does not occur in the dc-Sitter model, 
where one is free to choose the initial conditions at arbitrarily early times, like in the case of the Bunch-Davies state 
for which — > 0. 

We are particularly interested in the behavior of the quantum state at large times (fca/ao 3> 1). In this limit, Zk 
assumes the approximate form 



-2ik- 



2 ka 



1 + TZke 



-2ik- 



o 



(77) 



The above expression makes it clear that in the adiabatic late time limit, the mean particle number in any given mode 
saturates at a finite value: 



lim 

ka/ao — > 



[nk) 



1 



"0 



(78) 



The plot of (nk) (for the particular case where the instantaneous vacuum is chosen at ka/ao = 1/2) is shown in 
figure 7(a). Starting from zero, the mean number is characterized by transient oscillations that get steadily suppressed 
as the mode turns adiabatic. Figure 7(b) displays the corresponding complex trajectory of z^. which ends up going 
around in a circle centered on the origin of the complex plane, indicative of a finite, asymptotically constant particle 
number. These are identical in form to what wc found in the case of the constant electric field. 
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FIG. 7: Plots of the particle number (uk) and the complex trajectory of Zk as functions of the dimensionless variable ka/ao 
for the radiation-dominated model, for the case when the initial condition is imposed at ka/ao = 1/2. Sub- and super- Hubble 
regimes are color-coded blue and red respectively, (uk) saturates at a finite value at late times (sub-Hubble scales). Zk starts 
from zero, and at late times ends up circling around the origin of the complex plane. 



Let us look at the Wigner distribution in this case. It is again described in terms of the following two functions: 

l+7^;.e"'"^'^ 

^fc(a) 



and 



Jkia) 



(l-|7^fe|2)fca2 
2fca2lm (7^fce"^*''^) 



(79) 



1 + Tike 



-2ik- 



aoa 



(80) 
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and these are depicted in figures 8(a) and 8(b). The Wigner function clearly starts ofi^ in an uncorrelated state 
(J7/c = 0) at a = ai and with a finite spread of magnitude cr^ = fc~^a^^. However, at late times, erf ^ and this 
corresponds to a distribution sharply peaked along the p-axis. (This is presumably a consequence of the a — > oo limit 
of the scale factor.) As for the corresponding general classical phase space trajectory, it is given, in terms of the 
conformal time 77 by 



C 



cos{krj + c) 



V 



where C = Ce"^ is an arbitrary complex number. The corresponding momentum has the form 

Pk = = ^0,0^ [krism{kri + c) + cos{kri + c)] . 
Expressing qk in terms of pk yields the following expression for a general classical trajectory: 



Pk 



± 



Ck 



Pi 



alC^il + k^T]"^) 



(81) 



(82) 



(83) 



where, again, rj is to be treated as an implicit function of q^, making use of Eq. (81). Every trajectory at the initial 
moment starts from some point {qk{o,i),Pk{a.i)) in the phase space. But we are particularly interested in the late 
time behavior here. It is evident that in this limit, — > and we get a trajectory along the p-axis. This matches with 
the corresponding limiting form of the Wigner function. Clearly, we here have yet another instance of the Wigner 
function showing peaking on the classical trajectory despite the absence of continuing particle creation. 
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FIG. 8: Variation of the dimensionless functions a^ak/k and kj^/a^ for the radiation-dominated model. Sub and super Hubble 
regimes are color-coded blue and red respectively. The Wigner function starts uncorrelated at ka/ao = 1/2 with a finite value 
of al; at sub-Hubble scales, —* causing the Wigner function to get peaked on the p-axis. 



This may be compared with the late time variation of the classicality parameter, which is shown in figure 9. Sk 
is oscillatory with an amplitude less than unity; this of course is a consequence of the nature of the underlying q-p 
correlation, which, as can be directly figured using Eq. (79) and Eq. (80), remains finite and oscillatory over the course 
of the late time adiabatic evolution. Recalling the corresponding behavior in the non-adiabatic de Sitter case we dealt 
with earlier, the evolution, with regard to this feature, differs sharply in the two scenarios (as does the mean particle 
number) . 

The imaginary part of the time dependent effective action at late times directly expresses the asymptotic particle 
content in the quantum state, and an out-vacuum is definable in the adiabatic region. On the other hand, making 
use of the approximate expression for Zk in Eq. (77) gives the following expression for the real part of the effective 
Lagrangian at late times: 

m 

Except for the first term in the brackets, the remaining terms (up to first order in 1/a) are oscillatory in the a —^ 00 
limit, so their contribution to the effective action may be expected to be highly suppressed. Thus, the first non-trivial 
contribution to the effective action from the late time adiabatic evolution comes from the non-oscillatory term with 
magnitude aQ/{4ka'^), which is equivalent to e^{a)k/16. 
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FIG. 9: Plot of the classicality parameter Sk as a function of the dimensionless variable ka/ag for the radiation-dominated 
model, for the case when the initial condition is imposed at ka/ao = 1/2. Sub- and super-Ifubble regimes are color-coded blue 
and red respectively. Sk at late times ends up oscillating about zero with a finite [constant] amplitude. The variation of Sk 
reflects the complex trajectory of Zk, which at late times is a circle centered on the origin. 



Based on the analysis carried out above, combined with that of the de Sitter model, it is evident that peaking of 
the Wigner function on the corresponding classical phase trajectory at late times is a fairly general feature, seemingly 
independent of the adiabatic vs. non-adiabatic nature of the time evolution. In contrast, the classicality parameter, 
providing an alternative handle to quantify the variation of the Wigner distribution, shows appreciably dissimilar 
behavior in the two examples, and bears a rather close relation with particle production. As alluded to earlier, these 
observations are pretty much similar to the results obtained our in the toy models in Paper I revealed. 

The analysis of the preceding sections brings out the contrast of time evolution of quantized field modes in the 
inflationary scenario with that in the electric field case, and we would like to emphasize it once again here. A 
comparison of the plots for various quantities evaluated for the de Sitter case in section III B with those in section III A 
makes this point amply clear. In particular, a constant electric background provides adiabatic in as well as out regions 
where the notion of a particle can be unambiguously defined, which is not the case for the de Sitter model, where the 
late time evolution is highly non-adiabatic. Moreover, the power spectrum (related to the two-point function for the 
field) of field modes does not provide a direct measure of the particle content in the adiabatic hmit in the electric field 
scenario (see Eq. (32)) due to the asymptotically oscillatory nature of Zk; in contrast, for the de Sitter background, as 
is clear from Eq. (55) such a correspondence is made possible in the late time non-adiabatic phase of evolution where 
Zfc ^ 1. As pointed out before, the electric field case is qualitatively much closer to a decelerating universe like the 
radiation-dominated example, at least in being characterized by a late time adiabatic phase. The similarities in the 
nature of the plots in sections III A and IIIB make this fairly clear. 



4. A radiation-dominated universe going over to an asymptotically de Sitter phase 



We would like to provide another illustration of our formalism by considering quantized field evolution in a cos- 
mological model which decelerates (radiation-dominated) initially before making a transition to a de Sitter phase of 
accelerated expansion. The Friedmann equation for the scale factor describing such a universe [11] is given by 



fPrO , \ 



3 \ a 

(pro being the radiation energy density at a = oo = 1) and the solution in terms of the cosmic time t has the form 



(85) 



ait) 



Pro 
PA 



1/4 



sinh I 2\l PA t 



1/2 



PrO 
PA 
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(sinh(2T))'/' 



(86) 



where we have introduced the dimensionless variable r = \J^^-~PA t. 

The adiabaticity parameter corresponding to this expansion factor diverges in both extreme limits (a ^ and 
a —> oo), so asymptotic in and out states are not definable. Nonetheless, one can solve Eq. (4) to determine the time 
evolution of the quantum state of every mode. 
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FIG. 10: Plots of the particle number (nt), complex trajectory of Zk and classicality parameter Sk as functions of the di- 
mensionless variable r for the radiation+Lambda model. The instantaneous ground state is imposed at r = 0.001, and 
radiation-dominated and Lambda-dominated phases have been color-coded blue and red respectively, (nt) diverges at late 
times (super Hubble scales). Zk, after circling around for a while, asymptotically ends up at Zfc = 1. This is reflected in the 
evolution of Sk, which oscillates in the intermediate sub Hubble phase before saturating at unity. 



Figures 10(a), 10(b) and 10(c) show the numerical solutions for the mean number of particles, complex path of Zk 
and classicality parameter for a mode k corresponding to ^/S/SirGk/ {proPAo)^^^ — 10, plotted against the variable 
T. Here, the initial condition Zk — has been imposed at r = 0.001 (which corresponds to Rn/^p = 0.5) and the 
radiation-dominated and Lambda-dominated phases have been color-coded blue and red respectively. The average 
particle number, after a short period of transient oscillations (when the mode is temporarily inside the Hubble radius) 
monotonically increases at late times, and this mirrors the behavior of z^, which circles around the origin for a while 
before finally winding up at the point Zfc — 1. The behavior of Sk, saturating at unity after a phase of oscillatory 
variation, follows from this. In fact, one can see pretty clearly that the plots for this model can be thought of as 
obtained by 'patching up' the corresponding plots for the pure radiation-dominated and de Sitter examples which 
were considered earlier. Our time-dependent formalism thus provides a convenient way to visualize the differences in 
mode evolution during the decelerating and accelerating phases of expansion. 



5. Approximate analysis of accelerating models 



We would now like to verify that the features suggested by the preceding examples are fairly generic, i.e. they 
hold to a reasonable extent for all power law models (i.e. for all values of the index p). To this end, an approximate 
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analysis of the time evolution for arbitrary p is outlined below. We will first consider the range p > 1 corresponding 
to inflation, followed by the case of decelerating models with < p < 1. 
The accelerating models are described by the following scale factor [14] : 

a{t) = aofP {p > 1) (87) 

which, in terms of the conformal time coordinate can be re-expressed as 

a(r?) = {-nvy+'^ (-00 < 7? < 0) 



where 7 = {2p — 1)/(1 — p) with 7 < —2 and Ti, = [p — l)aj^^. De Sitter inflation is a particular example belonging to 
this set, and corresponds to the choice 7 = —2. Needless to say, all these models are characterized by non-increasing 
Hubble radius, and the evolution of every field mode turns strongly non-adiabatic at late times following Hubble 
crossing. 

We would like to consider a scenario where the quantum state of any mode starts out in the instantaneous vacuum 
in the asymptotic past. The solution to the equation for /Xfc, Eq. (5), which gives vanishing particle number in the 
7] ^ —00 limit is given by 

l,k{v) = N{kvrHi'\kr]) (89) 

where N = ^jT^jlk" and 7 has been replaced by the new variable v = —7 — 1/2; the normalization has been chosen 

so that /i'j*/ife — /i'j-Mfc = ^"i/o."^: so the Wronskian is set equal to unity. 

Fixing the form of pk, in principle completely specifies the time evolution of the state in all respects. Since we 
however are particularly interested in determining the late time behavior, the rj limit of the various quantities will 
be considered. Let us start with the solution for fj.). in Eq. (89). At super-Hubble scales {k\r]\ <^ 1), the approximate 
form of the Hankel function is given by 

Hi^\k,) . A('-^y\BC^y~^ (90) 

where 



2 V 2 



^ r(l - z/)sin7r!/ ' ^ r(2 - z/) siuTT!/ 

and so 

Pk{v)^N[2''A + 2-'-^B{krif + 0{rj^+^)] , fi'M ^ N [2''-'Bk^r] + 0{ri'+^)] (A > 0). (92) 

This gives the following approximation for the energy in the mode at late times (super-Hubble scales) to the lowest 
order of approximation: 



= ^W + kW) « 2'^-^nk\Amkvy-'' ^ ' J\ {kn)'-"- (93) 

^ (r(l — vj) sm TTU 



Since > 3/2 for this class of models (corresponding to p > 1), the energy in every mode diverges at sufiiciently late 
times {kri 0). This also means that the mean number of particles, being directly related to the energy, diverges as 
well: 

K Z K (1(1 — 1/)) sm TTJ/ 

Let us now move over to the phase space and consider the behavior of the corresponding Wigner function at late 
times; this is specified in terms of the functions and Jk, whose late time approximations have the forms 

al = 2\fik\^ « [22^|A|2 + 2^''-^{AB* + A*B){k7]f] (95) 

and 
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It is clear that cr^ has a finite limiting value, while J]^ is an increasing function of time, and diverges at super-Hubble 
scales. This represents a Wigner function that sharply peaks on the p-axis at late times. 

In contrast, in the asymptotic past as Zk — > 0, cr^ a — > oo while — > 0, describing a Wigner function that 
peaks around the g-axis. 

The behavior of the Wigner function needs to be compared with that of the corresponding classical trajectory in 
phase space. A general solution to the classical equation of motion is given by 



5fe = Re[Cr?''F«(fcr?) 



where C is an arbitrary complex constant. At sufficiently early times (j] 
Hankel function, it can be shown that 



(Ik 

while the conjugate momentum is given by 



Re 



C77 



(97) 

00), using the asymptotic form of the 



Pk = a^q'k 



Re 



k Im 



i i(fer,+0)\ 



(99) 



which implies that in the asymptotic past, pk ^ with an oscillatory qk whose amplitude sharply increases in the 
r] —00 limit, and the classical trajectory starts close to the q-axia. In contrast, at late times (A;|?7| <C 1), we have 



qk ~ R-e 

with the corresponding momentum being 




Pk = a^(/k 



2(-W)i-2'^Re CB 



n 



2-2v 



(100) 



(101) 



Since u > 3/2, qt/pk ^ as ^ 0, and we have a trajectory that ends up on the p-axis. Thus, based on the above 
limits one concludes that the Wigner function, for all p > 1, peaks on the classical trajectory in both asymptotic regions. 
This generalizes to all accelerating power law models what was earlier found for the de Sitter case, namely that this 
aspect (peaking on the classical trajectory) of the Wigner function may not provide much information regarding the 
approach to classicality. 

Furthermore, using the approximations for cr^ and J7fe in Eqs. (95) and (96), it trivially follows that the q-p c;orrelation 
(being just the product Jkcrl) will diverge as Jk ^ 00 at late times, implying a classicality parameter which ends up 
saturating at unity, accompanying strong particle production. This feature, too, was found to occur in the de Sitter 
example, and now has been generalized to all accelerating models of the type a cc . 

We will now move on to examine the other class, that of decelerating models, and explore the generality of the 
features the radiation-dominated model taken up earlier suggested. 



6. Approximate analysis of decelerating models 



These models are characterized by the scale factor 

a{t) = aotP (0 < p < 1) 
and in terms of conformal time, have the form 

a{n) = {-nr]y+^ (0 < r/ < 00) 



(102) 
(103) 



where 7 = (2p — 1) / (1 — p) and Ji = {p— 1)0^^. In contrast with the situation in accelerating models, the adiabatic 
parameter given by Ck = 2/kri is strongly suppressed in the late time sub-Hubble regime [krj 00). 

Like in the previous case, we would like to focus on the late time behavior of various quantities associated with the 
evolving quantum state. So we start with the general solution to Eq. (5), which is given by 



AkHi'\kr]) + BkHi^\kri) 



(104) 
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where v = —7 — (1/2). It should be recalled that the 77 ^ limit here is highly non-adiabatic. so a natural in- vacuum 
state is not definable in this region. We will therefore just work under the assumption that the state is chosen to 
coincide with the instantaneous vacuum at some moment t] = rji, which basically fixes B^/Ak- At sufficiently late 
times, however, every mode (in any arbitrary state) will eventually turn sub-Hubble and evolve adiabatically, so we 
have, for large values of krj, 



IJ'kiv) 



{AkB 



kf]- 



+ Bke 



ikr)-\-i 



/2fca(r?) 

which gives the following approximation for Zk at late times; 



) 



4:kr] 



(105) 



(106) 



where TZk = (6^/^^) exp (ztti/ + z7r/2). Using this approximate expression, one can derive approximations for the 
particle number as well as the correlation parameter in the adiabatic limit: 



(rik) 



1-W 



2kr] 



(107) 



which implies that the particle number saturates at the finite value IT^fcl^/ (l — I'^feP) = l-SfcP/ (|-4fc|^ — |Bfe|^), and 
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. (108) 



In the late time limit (fcry 00), only the first oscillatory term survives, and describes a correlation function that 
oscillates with a finite amplitude. Thus, as particle production is suppressed in the adiabatic limit, the q-p correlation 
remains bounded, generalizing the behavior which was found earlier in the radiation dominated model. 

From the expressions for a\ and Jk (see table), it is clear that in the strongly adiabatic limit [kr] — » 00) since \zk\ 
does not go to unity, the behavior of these functions is primarily determined by the late time limit of ruk = c?\ this 
implies cr^ ^ while Jk is oscillatory (due to the oscillatory nature of Zk) with an amplitude that is proportional to a^, 
describing a Wigner function that progressively gets peaked on the p-axis at late times. The Wigner function, however, 
was started off at 77 = 770 in an uncorrelated state [Sk = 0) that is fairly spread out (with al{r]i) = k~^a~^{r]i)). Thus 
the Wigner function gets sharply peaked with time even though the mean particle number remains finite all along. 
As for the classical trajectory, it has the same form as in the p > 1 models, Eq. (97); however, the late time limit here 
corresponds to taking — > 00, and this gives 



Re 



j(fer,+0) 



Ct]" 2 cos {krj + (l) + c) 



where C has been written as Ce"^, and 



Pk = a^q'k 



2v 



C 



T] ^ ^ cos [kr] + (f>+ c) — kr] sin [kr] + + c) 



(109) 



(110) 



Since < 1/2 for < p < 1, we have qk{ii) ^ as 77 ^ 00, and the classical trajectory (for arbitrary C) ends up 
on the p-axis, matching with the late time behavior of the Wigner distribution. As in the case of the classicality 
parameter above, this too is a generalization of what was seen in the radiation-dominated example. 

To conclude, based on the approximate analysis of the last couple of subsections, the two distinct late time limits of 
power law expansion, adiabatic; vs. non-adiabatic, show clearly differing behavior as regards particle production and 
variation of the classicality parameter (found to be sensitive to the particle creation) providing one possible measure 
of phase space correlations. Peaking of the Wigner function on the classical trajectory, on the other hand, is a fairly 
general feature, happening whenever one of the oscillator parameters blows up/goes to zero, and is quite independent 
of the particle production, i.e. adiabatic vs. non-adiabatic evolution. The observations drawn here, it may be noted, 
are broadly similar to, and thus to an extent corroborate, what we learned from the various toy models dealt with in 
Paper I. 



25 



IV. DISCUSSION 

In the previous sections, we have attempted to shed some Ught on a few conceptual issues arising in the study of 

the quantized time dependent oscillator, focusing our attention on a couple of conventional examples drawn from field 
theory in external backgrounds that hold particular physical relevance; this has been intended as an exploitation of 
the analysis (based on the Schrodinger picture formalism) that was carried out in Paper I. A quick recapitulation of 
the key features of the time evolution (of a given fourier mode) in the asymptotic limits, as revealed by our analysis 
of the examples of de Sitter and radiation models and the constant electric background case, is tabulated below: 



REGIME OF 
INTEREST 


dc Sitter model 
a{r,) = -1/Hrj 


Radiation dominated model 
a{r)) ^ aori 


Constant Electric field 
A^(t) - -Et 




Tj — ^ —OO 


ri ^ (But initial condition is 
chosen here at = 770 7^ 0) 


t —00 




Adiabatic 


Highly non-adiabatic 


Adiabatic 


EARLY TIMES 


{uk) 


{rik) =0 


(nk) 






(T^ : finite 






Ti. — ► 


7u — 


Ti, — * 






Sk = ^ 


5k ^ 




Strongly non-adiabatic 


Adiabatic 


Adiabatic 




?7 -> 


7} ^ 00 


t ^ 00 




(Uk) — > OO 


{nu) constant 


(nk) constant 


LATE TIMES 


(T^ — ^ finite 
constant 








Jk ^ OO 


J'k : oscillatory with 
increasing amplitude 


J7k : oscillatory with 
increasing amplitude 




Sk 1 


Sk '■ oscillates about Sk = 
with fixed amplitude 


^k : oscillates about Sk = 
with fixed amplitude 



For the examples (particularly those of the de Sitter and electric field models) considered here, the Wigner function 
for any field mode is found to get peaked on the corresponding classical phase space trajectory not just at late times, 
but even in the asymptotic past. This renders the interpretation of classicality based solely on peaking on the classical 
trajectory dubious, and naturally brings in the need to have an additional measure to specify the degree of correlation, 
if one seeks to address the question of approach to classicality confining oneself solely to the phase space evolution of 
the system. Our choice of the classicality parameter, which is found to closely track the particle production, proves 
fairly reasonable in this regard and is seen as having a more sensible interpretation in some situations, like being zero 
in a vacuum state even though the Wigner function may peak in this limit etc.. 

Let us reconsider the various limiting cases we came across in our examples in a step-by-step manner. First, let us 
consider the limit of early times. For all three cases considered, the function as well as the classicality parameter 
vanish at the moment when the mode is started off in the instantaneous ground state, and this clearly represents 
an uncorrelated Wigner distribution, although it is peaked in a limited region in phase space. (In the de Sitter and 
electric field cases, one can start off in the instantaneous vacuum state at t = — oo, but for the radiation dominated 
universe, a non-zero initial time rji is chosen.) 

As for the late-time behavior, the three examples differ in specific respects. For the de Sitter model, cr^ approaches 
a finite limiting value while Ji^ blows up. It can also be seen that the classicality parameter smoothly goes to unity 
in the a — ^ oc' limit [while it remains nearly zero at sub-Hubble scales in the asymptotic past]. Clearly, these features 
are indicative of a fairly correlated phase space distribution, and this is in line with the usual conclusion of classicality 
emerging on super-Hubble scales [10, 14, 17, 18]. For the radiation-dominated case, the Wigner distribution ends up 
getting peaked on the p-axis (coinciding with the corresponding classical trajectory) at late times on account of the 
fact that cr^ — * 0. 

As regards the variation of the classicality parameter, it turns out that this measure is zero to start with (at the 

instant when the mode is in the instantaneous ground state) but subsequently (for rj > tjq) develops an oscillatory 
profile that has a finite amphtude at late times. So, it can probably be argued that in some sense, the state gets 
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significantly more correlated in comparison with the initial state, although this is in contrast with the de-Sitter case 
[where 5k ^ 1 at late times], and is not so obvious. The difference is clearly tied with the adiabatic vs. non-adiabatic 
nature of the late time evolution. 

Finally, for the electric field example, the late time variation qualitatively proceeds along lines similar to the 
radiation-dominated model, in that cr^ 0, representing a Wigner function peaked on the p-axis [and which coincides 
with the classical trajectory]. However, the classicality parameter exhibits a variation that's more straightforward to 
interpret: it vanishes in the asymptotic past (<Sk ^ as i ^ — oo) but is oscillatory with a finite, nearly constant 
amplitude in the t oo limit. So in terms of iSk, the phase space correlations clearly appear to grow more significant 
in the far future in relation to early times. The de Sitter and electric field cases thus display a fairly unambiguous 
emergence of phase space correlations with time, but the evolution in the radiation model, in this particular respect, 
is a little less clear. The bottom line, thus, is that particle production is not necessary for peaking of the Wigner 
function on the classical trajectory as the three examples that were considered indicate, where this was found to 
happen even as (nk) saturated at late times or decreased to zero in the asymptotic past. (It may however be noted 
that one or the other of the oscillator parameters blows up/goes to zero in these limits.) Our analysis thus points 
toward the fact that peaking on the classical trajectory alone need not suffice to regard a system as turning classical; 
one needs to supplement this criterion with looking at some sensible measure of phase space correlations to settle this 
question. The classicality parameter that we have used is a simple yet useful construct in this regard, and leads one 
to conclude that strong particle creation, although not necessary for causing the Wigner function to get concentrated 
on the classical trajectory, is a necessary condition for driving the quantum state to become highly correlated (i.e. 
(qp) oo). 

ACKNOWLEDGEMENTS 
G.M. is supported by the Council of Scientific & Industrial Research, India. 



[1] V. Mukhanov and S. Winitzki, Inroduction to Quantum Effects in Gravity (Cambridge University Press, Cambridge, 

2007); N. D. Birrcll and P. C. W. Davics, Quantum Fields in Curved Space (Cambridge Univ. Press, Cambridge, 1982); 

S. A. Fulling, Aspects of Quantum Field Theory in Curved Spacetimc (Cambridge Univ. Press, Cambridge, 1989). 
[2] B. S. DeWitt, Pliys. Kept. 19, 295 (1975); D. N. Page, New ,J. Phys. 7, 203 (2005) [arXiv:hep-th/0409024]; T. Padmanab- 

han, Phys. Rept. 406, 49 (2005) [arXiv:gr-qc/0311036]; L. Sriramkumar, T. Padmanabhan, Int. J. Mod. Phys. D 11: 1 

(2002) [arXiv:gr-qc/9903054]; R. Brout, S. Massar, R. Parentani, P. Spindel, Phys. Rept. 260, 329 (1995). 
[3] J. Schwinger, Phys. Rev. 82, 664 (1951); for a text book discussion see C. Itzykson and J. B. Zuber, Quantum Field 

Theory (McGraw-Hill, New York, 1980). 
[4] T. Padmanabhan, Pramana, 37, 179 (1991); G. V. Dunne, in 'From Fields to Strings: Circumnavigating Theoretical 

Physics', Shifman, M. (ed.) et al. [arXiv:hep-th/0406216]; S. P. Kim et al., Phys. Rev. D 75, 045013 (2007); S. P. Kim, D. 

N. Page, Phys. Rev. D 73, 065020 (2006); Y. Kluger et al., Phys. Rev. Lett. 67, 2427 (1991); F. Cooper and G. C. Nayak, 

arXiv:hep-th/0611125. 

[5] C. Kiefer, Phys. Rev. D 45, 2044 (1992); C. Kiefer et al., Class. Quant. Grav. 8, L185 (1991); L. Sriramkumar and 

T. Padmanabhan, Phys. Rev. D 54, 7599 (1996). 
[6] L. Parker, Phys. Rev. Lett. 21, 562 (1968); Phys. Rev. 183, 1057 (1969); Phys. Rev. D 3, 346 (1971); L. P. Grishchuk, Sov. 

Phys. JETP 40, 409 (1975); B. L. Hu, G. Kang and A. Matacz, Int. J. Mod. Phys. A9, 991 (1994) [arXiv:gr-qc/9312014]; A. 

Campos and E. Verdaguer, Phys. Rev. D 49, 1861 (1994) [arXiv:gr-qc/9307027]; S. Habib, C. Molina-Paris and E. Mottola, 

Phys. Rev. D 61, 024010 (2000) [arXiv:gr-qc/9906120]. 
[7] Gaurang Mahajan and T. Padmanabhan, arXiv:0708.1233 [gr-qc]. 

[8] M. R. Brown, in 'Quantum Theory Of Gravity', Christensen, S. M. (ed.), 243-250 (1984); K. Srinivasan and T. Padman- 
abhan, Phys. Rev. D 60, 024007 (1999) [arXiv:gr-qc/9812028]. 

[9] R.H. Brandenberger, Rev. Mod. Phys. 57, 1 (1985); V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. 
Rept. 215, 203 (1992). 

[10] J. Martin, arXiv:0704.3540 [hep-th]; Lect. Notes Phys. 669, 199 (2005) [arXiv: hep-th/0406011]. 

[11] T. Padmanabhan, Theoretical Astrophysics, Volume III: Galaxies and Cosmology, (Cambridge University Press, Cam- 
bridge, England, 1999), Sec. 5.7. 

[12] A. R. Liddle and D. H. Lyth, Cosmological Inflation and Large-Scale Structure (Cambridge University Press, Cambridge, 
1999). 

[13] S. W. Hawking, Phys. Lett. B 115, 295 (1982); A. A. Starobinsky, Phys. Lett. B 117, 175 (1982); A. H. Guth and S.- 
Y. Pi, Phys. Rev. Lett. 49, 1110 (1982); J. M. Bardeen, P. J. Steinhardt and M. S. Turner, Phys. Rev. D 28, 679 (1983); 
A. D. Linde, Phys. Lett. B 116, 335 (1982); L. F. Abbott and M. B. Wise, Nucl. Phys. B 244, 541 (1984); T. Padmanabhan, 
Phys. Rev. Lett. 60 2229 (1988); T. Padmanabhan, T. R. Seshadri and T. P. Singh, Phys. Rev. D 39, 2100 (1989). 



27 



[14] L. Sriramkumar and T. Padmanabhan, Phys. Rev. D 71, 103512 (2005). 
[15] M. Morikawa, Phys. Rev. D 42, 2924 (1990). 

[16] T. Padmanabhan and T. R. Seshadri, Phys. Rev. D 34, 951 (1986); T. Padmanabhan, Class. Quant. Grav. 3, 911 (1986). 
[17] E. Wigner, Phys. Rev. 40, 749 (1932); A. L. Matacz, Class. Quantum Grav. 10, 509 (1993); L. P. Grishchuk and 

Y. V. Sidorov, Phys. Rev. D 42, 3413 (1990); T. P. Singh and T. Padmanabhan, Annals of Physics 196, 296 (1989); 

D. Polarski and A. Starobinsky, Class. Quantum Grav. 13, 377 (1996) [arXiv:gr-qc/9504030]; S. Habib and R. Laflamme, 

Phys. Rev. D 42, 4056 (1990). 

[18] E. Calzetta and B. L. Hu, Phys. Rev. D 52, 6770 (1995) [arXiv:gr-qc/9505046]; D. H. Lyth and D. Seery, arXiv:astro- 
ph/0607647v2; C. Kiefer and D. Polarski, Annalcn Phys. 7, 137 (1998) [arXiv:gr-qc/9805014v2]; M. Mijic, Phys. Rev. D 
57, 2138 (1998) [arXiv:gr-qc/9801094vl]; J. H. Heyl, arXiv:gr-qc/0602079vl. 

[19] I. S. Gradshteyn and I. M. Ryzhik, TabJe of Integrals, Series and Products (Academic Press, London, 1994). 

[20] G. B. Arfken and H. J. Weber, Mathematicai Methods for Physicists (Academic Press India, 2001). 



